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ABSTRACT 

We have analyzed the observed image separations of the gravitationally lensed 
images of QSOs for a possible correlation with the source redshift. Contrary to 
the previously noted anti-correlation based on a smaller data set, no correlation 
is found for the currently available data. We have calculated the average image 
separations of the lensed QSOs as a function of source redshifts, for isothermal 
spheres with cores in a flat universe, taking into account the amplification bias 
caused by lensing. The shape of the distribution of average image separation as 
a function of redshift is very robust and is insensitive to most model parameters. 
Observations are found to be roughly consistent with the theoretical results for 
models which assume the lens distribution to be (i) Schechter luminosity function 
which, however, can not produce images with large separation and (ii) the mass 
condensations in a cold dark matter universe, as given by the Press-Schechter 
theory if an upper limit of l-7xl0 13 M© is assumed on the mass of the conden- 
sations. 

Subject headings: cosmology: gravitational lensing 

1. Introduction 

The phenomenon of gravitational lensing is extremely useful for understanding the 
large scale structure of the universe. Studying the observed properties of the images in 
individual lenses can provide knowledge of the mass distribution and mass to light ratio 
in these lenses (Grogin and Narayan, 1996; Tyson, Kochaski and Dell'Antonio, 1996). A 
statistical analysis of the observed lens systems on the other hand can help us to restrict 
the values of the cosmological parameters (Falco, Kochanek and Munoz, 1998; Link and 
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Pierce, 1998). Based on the frequency of gravitational lensing, upper limits have been 
placed on the value of the cosmological constant (Fukugita et al 1992, Kochanek, 1996). 
It has recently been pointed out by Park and Gott (1997, hereafter PG) that the image 
separations of the observed lenses are strongly negatively correlated with source redshifts, z s , 
the anti-correlation being much stronger than that predicted by standard cosmologies. PG 
considered several possible ways to strengthen the anti-correlation. A steeper mass profile, 
merger of galaxies and increase in their mass by in-fall strengthen the anti-correlation, 
but not sufficiently to explain the observations. PG considered point masses and singular 
isothermal spheres (SISs) for the lenses and assumed all the lenses to be identical. They 
also, did not consider the effect of amplification bias in their calculation of average image 
separation. 

Williams (1997, hereafter LLRW) showed that the theoretical upper limit to the image 
separations shows a strong anti-correlation with the source redshift, consistent with the 
observations, provided (i) the lensing galaxies have logarithmic surface mass densities that 
gradually change with radius, (ii) there is a dispersion in the lensing properties of galaxies 
like the central surface mass density or velocity dispersion and (iii) the characteristic length 
scale of dark matter halos of galaxies scales as L a , with a~ 0.4. 

More data on lensed QSOs have become available since 1997 and the number of lensed 
QSOs is now roughly twice the number in the PG sample. The PG sample also contained 
some doubtful lenses. The new data set has to be examined afresh for the presence or 
absence of the correlation between image separation and source redshift. The main aim of 
the present paper is to analyze the currently available data to look for a possible correlation. 
We then want to compare the observed distribution of average image separation as a 
function of source redshift with the results of a detailed calculation of statistical lensing of 
galaxies, with a view of obtaining constraints on the values of various parameters entering 
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the calculations. The plan of the paper is as follows. In section 2 we present the analysis 
of the data using various statistical tests to quantify the presence of any correlation. 
Calculations for lensing statistics in particular of the average image separation as a function 
of source redshift, for flat universe, for realistic lens distributions, taking into account the 
amplification bias, for different cosmological models are presented in section 3. In section 4 
and 5 we present the results and conclusions respectively. 

2. Is there a correlation? 

The anti-correlation noted by PG is, as observed by them, mainly due to five low 
redshift (z s <2.15), large separation (< A9 > greater than 4") images, namely, 0240-343, 
0957+561, 1120+019, 1429-008 and 2345+007. Out of these, only one (0957+561) was a 
confirmed lens. Since then Kochanek, Falco and Munoz (1999) have argued against the 
lensed nature of most of the wide separation QSO pairs (hereafter WSQPs). Studying the 
optical and radio properties of these pairs they concluded that all of the WSQPs, with 
< A9 > between 3 " and 10 " are binary quasars rather than gravitationally lensed images 
with a one sided 2 a (1 a) upper limit of 22% (8%) on the lens fraction of these QSO 
pairs. They explained the high incidence of occurrence of binary QSOs (which is 2 orders 
of magnitude higher than that given by the quasar-quasar correlation function) as being 
due to the enhanced quasar activity during the merger of galaxies. No lensing galaxy could 
be found in front of these WSQPs so that if they are lensed images, the lensing mass, 
though of the magnitude of a cluster, has to be completely invisible. Peng et al (1999) have 
searched for a lensing galaxy in front of Q 1634+267 at the lens optical wavelength using 
NICMOS and showed that the lens has to have M/L > 690 h 65 (1200 h 65 ) for Q =0.1 (1.0) 
and H = 65 h 65 km s" 1 Mpc -1 . They therefore, suggest that the double "images" may be 
binary QSOs rather than multiple images of a single QSO. Very similar spectra of the two 
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"images", however, defies this conclusion. Peng et al (1999) have compared the spectral 
similarities of 14 pairs of QSOs with separations between 3 " and 10 " which are having 
very similar redshifts (with the velocity differences between pairs being < 500 km s _1 ) with 
the spectral similarities of randomly chosen QSO pairs from the LBQS. They conclude that 
there is < 3% probability that a randomly drawn sample of 14 QSO pairs would show as 
similar spectra as the observed pairs. So unless a viable theory to explain the similarities 
of the spectra of QSOs in merging galaxies can be developed, one can not discard the lens 
hypothesis for the WSQPs. 

New data on lensed QSOs have become available since 1997. 49 confirmed or likely 
lenses have been compiled by the CfA/Arizona Space Telescope Lens Survey (CASTLES) 1 , 
out of which the source redshift is available for 39 QSOs. These include 15 lenses from the 
PG sample. For reasons stated above we may combine these with the 5 WSQPs from the 
PG sample (0240, 1120, 1429, 1634 and 2345) to get an extended sample of 44 lenses. We 
have considered 4 separate samples for our analysis, namely (1) PG20 ; the sample used by 
PG of 20 QSOs (2) PG15 ; PG sample without the 5 WSQPs (3) CAST39; the CASTLES 
sample of 39 lenses and (4) EXT44; the extended sample containing 39 CASTLES lenses 
and 5 WSQPs. The angular separation vs source redshift for these samples are plotted 
in Figures la, lb, lc, Id respectively. These samples were searched for the presence of 
correlation by performing Spearman rank correlation test. KS test was also performed to 
determine if the image separations for sources with redshift smaller than 2.5 and sources 
with redshift greater than 2.5 are taken from the same distribution. We also obtained 
best straight line fits to the unbinned data of the four samples which are plotted in the 
figure. The results are shown in Table I. It is clear that PG20 shows highly significant 
anti-correlation between < A9 > and z s . The anti-correlation weakens but persists at about 



1 see http:/ /cfa- www.harvard.edu/castles 
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1.8 a level after the removal of the 5 WSQPs (PG15). CASTLES data on the other hand 
do not show any correlation. Adding the 5 WSQPs to CASTLES data (EXT44) gives rise 
to anti-correlation but it is very weak and is statistically insignificant. KS test also shows 
that for EXT44 sample the probability that the lens separations for sources with redshift 
smaller than and larger than 2.5 are taken from the same distribution is quite large. We 
thus conclude that the strong anti-correlation noted by PG was mainly due to the inclusion 
of 5 questionable QSOs which formed a large fraction of about 25% of their sample. The 
present data of 39 confirmed lenses does not show any correlation by itself and even after 
the inclusion of the 5 questionable lenses which now form only 11% of the sample. 



3. Details of calculations 

As seen above no correlation exists between the image separations and the source 
redshifts and it seems possible that results of the standard theoretical models may be 
consistent with the observations without the need of any drastic assumptions. It is possible 
that we may be able to constrain some of the parameters entering the calculations. In 
this section we present the details of our calculation of lensing statistics which differ from 
previous calculations only in the way the amplification bias has been taken into account. 

We have assumed the mass distribution of the lenses to be that of ISCs as found to 
be required by Maoz and Rix (1993) in order to explain the lensing events in the HST 
snapshot survey. We have not taken into account the effects of ellipticity in the lensing 
galaxies. As the ellipticity mainly affects the relative numbers of double and quadruple 
lenses (Keeton et al, 1997) it will not affect our results. For the lens distribution as a 
function of mass and redshift we consider the Schechter luminosity function of galaxies 
and the mass condensations in a cold dark matter universe obtained using Press-Schechter 
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theory (hereafter referred to as the PS distribution). The former is given by 

<$>{L)dL oc (L/L*) a e- L/L *dL, (1) 

$ (L) being the number of galaxies with luminosity L per unit comoving volume. Here, a 
=-1.1 (Marzke et al, 1998). We have assumed the comoving number density of galaxies 
to be independent of redshift. The circular velocity dependence of the luminosity is taken 
to be L/L* = (a/a*) 4 , <7* = 225 km s" 1 for elliptical galaxies (Faber & Jackson, 1976, de 
Vaucouleurs & Olson, 1982, Kochanek, 1994) and L/L* = (a/a*) 2 - 6 , a* = 144 km s" 1 for 
spiral galaxies (Tully & fisher, 1977, Fukugita and Turner, 1991), which are assumed to be 
70 % of all the galaxies (Postman and Geller, 1984). The PS distribution is given by 

_ -3(1.67)^(1 + Z p dtoA -3(1 + ,)' 

n(v c ,z)dv c - (27r)3 /2 v 4 A(ro) (ll„ ¥t X expl 2A s (rJ ' " [ > 

Here n(v c , z) is the number density of mass condensations with circular velocity v c at 
redshift z, 5 C = 1.68 and the functional form of A(r Q ) for the CDM power spectrum of 
density perturbation is 

A(r D ) = 16.3 b" 1 (1 - 0.3909 r"- 1 + 0.4814 r°- 2 )~ 10 , (3) 

where b is the bias parameter. The mass and circular velocity of a halo are related to the 
comoving radius r D and redshift z by, 

M = ^ Po r 3 , V c = 1.67 (1 + zf 2 H r , (4) 

p being the mean density of the universe. The density fluctuation amplitudes are taken 
from N-body simulation work (Narayan and White, 1988, Mo,Miralda-Escude and Rees, 
1993). 

We consider several flat world models with different values of the cosmological constant. 
Following Fukugita et al (1992), we use the angular diameter distances between the lens 
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and the observer (Dol) between the source and the observer (Dos) and between the lens 
and the source (D LS ) and take the critical impact parameter for SIS lens to be 

a cr = 4tt(^) — , (5) 

c u os 

a\\ being the one component velocity dispersion for the lens equal to v c /V2. The angular 
diameter distance formulae for various types of world models are taken from Fukugita et al 
(1992). The lensing cross-section for ISC is given by (Hinshaw and Krauss 1987) a = ttIq, l 
being the maximum impact parameter for lensing, given by 

l = \{a 2 CT + 5a cr r c - 0.5r c 2 ) - Q.hrl' 2 {r c + 4a cr ) 3/2 ] 1/2 , (6) 

r c being the core radius of the lens mass distribution. The amplification of an image is 
obtained from A = ~. Here / is the impact parameter and b is the image position in the 
lens plane obtained by solving the lens equation 

b 3 + lib 2 + b{l 2 + 2a cr r c - a 2 cr ) + 2la cr r c = 0. (7) 
The average image separation for a given value, z s , is given by 

Z s OO lo 

j dzi j dv c J dl n s {z s , z b v c , M% m )A6(z s , z h v c ,l) 

< A9 >= ° 2 °- . (8) 

/ dzi J dv c J dl n s (z s , z h v c , M l ^ m ) 



Here, A9(z s , zi, v c , I) is the separation between the two brighter images produced by a lens 
at z; with circular velocity v c and impact parameter I and n s (z s , zi, v c , Mg m ) is the number 
of observable sources (with m b s <mu m ) at the redshift z s which will be lensed by galaxies 
at a redshift z\ with circular velocity v c and impact parameter /, such that both of the 
brighter images will have luminosity higher than that for M^ m . This can be obtained 
multiplying the QSO luminosity function by the optical depth for lensing and integrating 
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over the observable magnitude interval. This can be written as 

M l J l m +2.5logA 

n s (z s ,z h v c ,M< B im ) = J dM B <P{M B )^^^. (9) 

M B =— oo 

Here, M 1 ™ is the absolute magnitude of a source at z s corresponding to the limiting 
apparent magnitude m l B m of the survey, assumed here to be 18. The results are quite 
insensitive to the value of m B m . A(z s , zi, v c , I) is the amplification of the weaker of the two 
bright images and 0(M#) is the QSO luminosity function, which is taken from Wallington 
and Narayan (1993). 9 ^qI'qi^ is the optical depth for lensing of a source at z s by lenses 
with velocity dispersion v c at redshift z\ with impact parameter /. This is given by 

d 3 r(z s ,zi,v c ,l) cdt 

d Zl dv c dl =Mvc 4)2*1-, (10) 

n ; (t> c , z\) being the number density of the lenses with circular velocity v c at z; and 2 n I dl 
being the crosssection for lensing for impact parameters between / and dl. 

The unnormalized probability of lensing for a given image separation A# can be 
computed from 

z„ In 

d P (z s ,A 6) 
dA9 

where, 



= J dz x J dlm a (z a , Zl ,M l * m ,l,Ae), (11) 



M l ^ m +2.5logA 

m.(z.,z l ,Mlr,l,M)= J ^M g 0(M B ) g3r ^^f g) ( 12 ) 

M B =-oo 

is the observable number of QSOs (obtained as above by multiplying the QSO luminosity 
function with the optical depth of lensing and integrating over the observable magnitude 
range) lensed by lenses at redshift z t with impact parameter / producing image separation 
A9, A(z s ,z/, /, AO) is the amplification of the weaker of the two bright images with image 
separation A# and 

d 3 r(z s , zi,l, AO) dv c dt 

dz ld A6dl = ni{Vc > Zl) ^ l dA0 C dz~C (13) 
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v c being the circular velocity of lenses at z\ which will yield the image separation A8 
for impact parameter /. For this calculation we have ignored the dependence of image 
separation on impact parameter which is very weak (Hinshaw & Krauss, 1987) and have 
used the value of image separation by the ISCs at zero impact parameter. We have verified 
that this assumption does not lead to errors larger than 1 %. This is given by 



Ad = 2-^(1.0 (14) 
Dql a c ' 



>'cr 



The value of a cr for given A9 obtained from the above equation is used to obtain the 
necessary value of v c as 

vc = c L*Z D °s ) 1/2 - (15) 
lixD OL D LS 



4. Results and Discussion 

We have plotted in Figure 2a the unnormalized differential probability ^rj^ as a 
function of image separation for the two lens distribution functions for several redshifts. 
We have used a value of 0.2 kpc for the core radius for these calculations. This is the 
upper limit obtained by Wallington and Narayan (1993) from the observed absence of 
central images for the lensed QSOs. Models of individual lenses imply somewhat lower core 
radii (Kochanek, 1995). An upper limit of 1.4 kpc has been obtained from the results of 
N body simulations of gravitational collapse of density peaks by Dubinski and Carlberg 
(1991). The slope of the lens distribution (as a function of luminosity or circular velocity) 
for Schechter luminosity function is steeper than that for the PS distribution. As a result 
the probability distribution is broader for the PS distribution compared to that for the 
Schechter distribution. The slope of the probability distribution, for each lens distribution 
is, however, almost independent of the source redshift for A9 > 0.4" which indicates that 
the average value of A9 may not be very sensitive to z s . 
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Figures 2b and 2c show the effect of varying the values of various parameters on 
the differential probability distribution. Figure 2b shows that increasing the core radius 
suppresses the probability for small values of separations (curves 2,3,4,5). However, again, 
the slope of the distribution is independent of the value of the core radius for AO > 0.3" 
for r c < 0.2 kpc. The image separation increases with the increase in core radius as the 
probability for small separations is suppressed (Hinshaw and Krauss, 1987). However, the 
effect of including the amplification bias (which increases the probability for small values of 
A6 as can be seen by comparing curves 3 and 6) cancels the effect of the increase in the core 
radius as has been noted by Hinshaw and Krauss (1887). This results in the probability 
distribution being almost independent of r c for AO > 0.3" for the assumed range of r c 
values. Higher values of core radius change the probability distribution more significantly 
for larger values of separations. The effect of increasing the bias parameter is to reduce the 
probability for large values of AO as seen by comparing curve 1 and 3. In Figure 2c we have 
plotted the probability distribution for three cosmological models, including the cosmic 
concordance model proposed by Ostriker and Steinhaardt (1995). For PS distribution the 
peak shifts to lower A6 values with increase in A. 

The average image separation as a function of redshift has been plotted in Figure 3 
for several values of parameters. In this figure we have also plotted the observed image 
separations and the best fit straight lines for the EXT44 sample (small dashed line) and 
for the CAST39 sample (long dashed line). It is well known that the PS distribution tends 
to over predict large separation lenses (Kochanek, 1994; Flores & Primack, 1996). As a 
result the predicted average image separations are large compared to the observed values. 
In plotting the results in Figure 3 we have used a cutoff on the mass of the condensations. 
The probability distribution for a cutoff of 600 km s" 1 is shown in Fig2a for comparison. 
As is expected, the absence of large masses reduces the probability for large separations 
drastically. Porciani et al (2000) have shown that the probability of lensing is consistent 
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with observations if one assumes that the condensations with mass > 3.5 x 1O 13 M0 have 
non-singular mass distribution. However, as explained above, for r c < 0.2 kpc the effect of 
finite core radius is partly compensated by the amplification bias for < A9 » 0.3" and it 
is not possible to suppress the probability for separations > 4". A cutoff on the mass of 
condensation is therefore needed to reduce the probability for these separations. 

Figure 3a shows that the < A9 > values are almost independent of the source redshift. 
The decrease of < A9 > with z s is somewhat more pronounced for the assumption of PS 
distribution for small values (< 2) of z s . The decrease is larger when the distance formulae 
for empty beam are used compared to that for the case of filled beam (curves 1 and 2). 
The values of < A9 > are sensitive to the values of the upper limit on the circular velocity 
of the condensations used in the calculation. An increase in the upper limit from 750 km 
s _1 to 1000 km s" 1 increases the values of < A9 > by a factor >1.5 (curves 2 and 4). 
Change in the bias parameter (curves 4 and 6) again changes the absolute values but does 
not alter the slope of the distribution. As expected from Fig 2b, the values of < A9 > are 
almost independent of the core radius (curves 4 and 5). We thus see that the observed 
image separation as a function of redshift (as seen from the best fit line) is in reasonable 
agreement (in view of the uncertainties of the best fit given in Table I) with the theoretical 
results for an upper limit on circular velocity of about 600-750 km s" 1 for b=l for PS 
distribution. This value of course is not absolute as higher values will be needed for higher 
bias parameters. Use of Schechter luminosity function, however, produces considerably 
lower values of < A9 >. 

In Figure 3b we have shown the results for different assumptions about the change in 
the core radius with redshift and with circular velocities of the condensations as suggested 
by LLRW. Assuming r c oc -yc 5 ' 10 ' 2 (curves 4, 5 and 6) does not change the results 
significantly. Similarly the assumption of a redshift dependence of r c oc (1.0 + zi)^ 1 ^ 2 
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(curves 1 and 2) also changes the results only by a small amount. This is expected in view 
of the very weak dependence of < A6 > on core radius. Thus the observed absence of 
correlation is consistent with the results for a flat universe even if the core radius varies 
with luminosity contrary to the suggestion of LLRW. In Figure 3c we have plotted the 
results for different cosmologies. Values of < A9 > decrease with increase in A. 

The results presented above show that the average image separation is not very 
sensitive to (i) the values of core radii (ii) different assumptions about the dependences of 
core radii on luminosity or redshift and (iii) the values of cosmological constant between 
and ~0.65. The values are sensitive to the upper limit used for the circular velocity for 
the PS distribution. For b= 1, core radii between and 0.2 kpc and cosmological constant 
between and ~ 0.65, an upper limit of about 600-750 km s _1 on the circular velocities of 
the dark matter halos for the PS distribution is roughly consistent with the observations 
of image separations. It may be noted that the observed values of image separations 
do show considerable scatter at any given source redshift. A large scatter is expected 
from the flat probability distributions as a function of A6 (Figure 2). The differential 
probability of lensing as a function of angular separation can not be directly compared 
with the observations due to the small number of observed lenses at any given redshift. 
However, if we take all the observed lenses together, irrespective of their redshifts, then we 
see that 6 out of 49 lenses in the full CASTLES data and 11 out of 54 in the extended data 
including the WSQPs have image separation between 3 " and 8 ". This requires the ratio of 
probabilities for separations between 3 " and 8 " to that for separations between 0.3" and 3 
" to be 0.14 and 0.26 respectively for the two samples. We have calculated this probability 
ratio to be 0.62, 0.35 and 0.15 for upper limits of 750, 600 and 500 km s _1 respectively 
for the PS distribution for b=l, f2=l and r c =0.2 kpc. The ratio is 0.04 for the Schechter 
function. 
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Keeton, Christlein and Zabludoff (2000) have considered the detailed galactic luminosity 
function dependent upon type and environment (Bromley 1998a,b). They obtained a good 
match between their calculated probability and the image separation distribution for 49 
QSO lenses in the CASTLES data. We have calculated the probability distribution and 
the average image separation using the luminosity function used by them. The results are 
shown in Figs lc and 3c. The probability peaks at a higher value of image separation and as 
a result the values of < A9 > are higher and can be considered to be in agreement with the 
CAST39 sample. However we note that even for this luminosity function it is not possible 
to get values of image separations larger than 6 ". The ratio of probabilities for separations 
between 3 " and 8 " to that for separations between 0.3" and 3 " is 0.12. 

5. Conclusions 

We have analyzed the observed image separations of lensed QSOs for a possible 
correlation with the source redshift. A correlation was earlier noted for the then available 
data by PC The present data of 39 confirmed or likely lensed QSOs, even when combined 
with 5 wide separation QSO pairs which are doubtful lenses, do not show any statistically 
significant correlation. LLRW had shown that if the core radii scale with luminosity as L a4 
then an anti-correlation is expected theoretically because of the presence of dispersion in 
lensing properties of galaxies . We have calculated the average image separations of the 
lensed QSOs as a function of source redshifts, for isothermal spheres with cores in a flat 
universe, taking into account the amplification bias caused by lensing. We do not find 
the strong anti-correlation stipulated by LLRW. In fact the shape of the distribution of 
average image separation as a function of source redshift is very robust and is insensitive 
to the change in parameters. As a result we are unable to obtain meaningful constraints on 
any of the parameters. The assumption of the Schechter luminosity function is unable to 
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produce separations larger than 6 ". The use of PS distribution on the other hand yields 
large number of wide separation lenses even for non-zero core radius and nessecitates the 
assumption of a cutoff on the mass of the condensations or large values of core radii for 
condensations with large masses. It may be noted that the PS distribution has also been 
found to be remarkably successful in explaining the observed distributions of the QSO 
absorption lines (Das & Khare, 1999). 

The author is deeply grateful to the referee for his/her extensive, in-depth and 
educative comments which have brought the paper to the present form. 
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Fig. 1. — Image separation vs source redshift for four samples: (a) PG data of 20 QSOs (b) 
PG data without the 5 WSQPs (c) CASTLES data of 39 QSOs and (d) CASTLES data 
with the 5 WSQPs. Triangles represent the WSQPs. The best fit straight line for each data 
set is plotted. 

Fig. 2.- 

Fig 2a: Unnormalized differential probability of lensing dl ^ Ag ^ as a function of angular sepa- 
ration of the bright images. Solid and dashed lines are for the assumption of PS distribution 
and Schechter luminosity function respectively for the lens distribution. Curves labeled 
1,2,3,4 and 5 are for source redshifts of 1, 2, 3, 4 and 5 respectively. Dash-dotted line is for 
PS distribution, for z s =3 assuming a cutoff of 600 km s^ 1 on the circular velocities. All 
curves are for r c =0.2 kpc, f2=1.0, b=l and filled beam. 

Fig 2b: Unnormalized differential probability of lensing dl ^ Ae ^ as a function of angular sep- 
aration of the bright images, for several values of parameters. Lines labeled 2,3,4,5 are for 
r c =0, 0.1, 0.2 and 0.5 kpc respectively. Curve 6 has ignored the amplification bias. All these 
curves are for b=l. Curve 1 is for b=2 and r c =0.2 kpc. All curves are for Q—l, filled beam 
and PS distribution. 

Fig 2c: Unnormalized differential probability of lensing dl ^ A0 ^ as a function of angular sepa- 
ration of the bright images. Solid and long dashed lines are for PS distribution and Schechter 
luminosity function for the lens distribution respectively. Lines labeled 1,2 and 3 are for A= 
0, 0.65 and 0.9 respectively. Short dashed line is for the luminosity function used by Keeton 
et al(2000) for A =0.65. All lines, except for the short dashed line, are for r c =0.2 kpc, b=l 
and filled beam; short dashed line is for r c =0. 

Fig. 3.- 

Fig 3a: Average image separation as a function of source redshift. Curve 8 is for Schechter 
luminosity function for empty beam. All other curves are for PS distribution. Curves 1 and 
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2 are for filled and empty beams respectively assuming an upper limit of 1000 km s -1 on 
circular velocity. Curve 3 is obtained by ignoring the amplification bias. Curve 5 is for SIS 
lenses. All except curve 6 are for a bias parameter of 1. Curve 6 is for a bias parameter of 
2. All the curves are for £1=1.0. Curves 3-6 are for an upper limit on circular velocity of 750 
km s _1 while curve 7 is for an upper limit of 600 km s^ 1 . Curves 3-8 are for empty beam and 
all except curve 5 have r c =0.2 kpc, curve 5 is for r c =0. The observed image separations and 
the best fit straight line (dashed line) for these are also plotted. The WSQPs are represented 
by triangles. 

Fig 3b: Average image separation as a function of source redshift for PS distribution for 
different assumptions regarding the core radius. Curve 1 is for r c oc (1 + zi)~ 2 ; curve 2 is for 
r c oc (1 + Zi) -1 ] curve 3 is for constant r c ; curve 4 is for r c oc v °' 5 ; curve 5 is for for r c oc v c ; 
and curve 6 is for r c oc v\. All curves are for f2=1.0, r c =0.2 kpc (for zi=0 for curves 1 and 
2 and for v c =200 for curves 4,5 and 6), b=l, upper limit on circular velocity of 750 km s _1 
and empty beam. The observed image separations and the best fit straight line (dashed line) 
for these are also plotted. The WSQPs are represented by triangles. 

Fig 3c: Average image separation as a function of source redshift. Curves labeled 1,2 and 3 
are for A= 0, 0.65 and 0.9 respectively for PS distribution assuming an upper limit of 750 km 
s _1 on the circular velocity. Lines marked KE and KF are for the luminosity function used 
by Keeton et al (2000) for empty and filled beams respectively for r c =0. All other lines are 
for r c =0.2 kpc and empty beam. All lines are for b=l. The observed image separations and 
the best fit straight line (dashed line) for these are also plotted. The WSQPs are represented 
by triangles. 



-20- 



Table 1. Results of statistical tests and best fit 



Sample 


Spearman rank correlation test 


KS test probability 


Best fit parameters 


P 


a 


p/a 


a 


b 


PG20 


-0.57 


0.23 


-2.49 


1.87(-2) 


6.07±1.2 


-1.29±0.47 


PGf5 


-0.47 


0.27 


-1.78 


0.14 


3.84± 1.02 


-0.68± 0.37 


CAST39 


6.52(-2) 


0.16 


0.40 


0.53 


2.25± 0.60 


-0.13± 0.24 


EXT44 


-2.56(-2) 


0.15 


-0.17 


0.39 


3.23± 0.75 


-0.38± 0.31 
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